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ABSTRACT 

We perform a number of numerical simulations of the solar corona with the aim to understand 
how it responds to different conditions in the photosphere. By changing parameters which govern the 
motion of the plasma at the photosphere we study the behavior of the corona, in particular, the effects 
on the current density generated. An MHD code is used to run simulations, using a 20 x 20 x 20 Mm 3 
box with time spans ranging from one hundred to several hundreds of minutes. All the experiments 
show a fast initial increase of the current density, followed by a stabilization around an asymptotic 
value which depends on the photospheric conditions. These asymptotic average current densities as 
well as the turn-over points are discussed. 

Subject headings: sun: corona — sun: magnetic fields — stars: magnetic fields — stars: coronac 
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1. INTRODUCTION 

Over the last decades the solar corona has been, 
and continues to be scrutinized and studied in a 
search for th e expl a nation of it s heat i ng mechanism 
dRosner et al.1 Il978t iGolub et all 119801 : iHendrix et all 



119961: lAschwanden et all |2001[ ). With the advent of 



the computer, simulations become feasible and the re- 
production of observations through models have ad- 
yanced our understanding of the 1 MK solar corona 
(|Gudiksen fc Nordlundl l2005bh . A mong the theoretical 
work, two main ideas have become popular to explain the 
hot corona. One supporting a wave heating mechanism 
introduced by lAlfve n (1947j) where the energy is carried 
through magncto-hydrodynamic waves; and the other 
based on electric currents developed by th e movements o f 
the magnetic field lines, first suggested bv lParke"rl (|1972f ). 

It has been clear from the beginning that the mag- 
netic field play a very important role in the process of 
coronal heating, besides passively shaping coronal struc- 
tures. It works as an active agent carrying energy from 
the convective outer layer of the Sun, and participates in 
the energy storage and dissipation mechanism. The con- 
tinuous convective motions shuffle the footpoints of the 
magnetic field lines, twisting the magnetic field lines and 
increasing the complexity of the magnetic field configu- 
ration. In the context of the Parker nanoflare model, the 
kinetic energy is transmitted to the fluid and converted 
into heat by a dynamical non-equilibrium and, in spite 
of the very high electrical conductivity of the plasma, 
the corona is heated through the dissipation of current 
sheets. Several studies and observations have been made 
supporting these ideas, associating the coronal heating 
with the stress of magnetic field by the surface turbu- 
lence and showing a direct relation between the coronal 
emission and the location of the magnetic fields on the 
surface (Ivan Spevbroeck et al.lll970l IVaiana et al.lll973l : 
iKrieger et al.lll971lll976l to mention examples). 

It is known that main-sequence stars with masses 
AI < 1.3 M© (spectral types F and later) have a nu- 
clear burning core, surrounded by a shell with radia- 
tive transport and an outer convective envelope where 
the energy is carried out of the star by convective mo- 



tions. The convection zone is narrow in extent in F-type 
stars, but in the Sun (G2 V) it has a depth of 0.287 so- 
lar radii (jChristensen-Dalsgaard et al.lll99l[ ) and, with 
decreasing effective temperature the convective zone oc- 
cupies a larger and larger fraction of the stellar volume 
with late M stars being fully convective. The photo- 
sphere lies above this convective layer and marks a tran- 
sition between a medium where energy is carried by con- 
vection to a radiative efficient domain. The granular 
structure at the stellar surface is fundamentally different 
from the flow topology below the photosphere. Numer- 
ical models suggest that convection is mostly driven by 
the downflows due to ra diative cooling near the surface 
(|Stein fc Nordlundl 1 19981 ) and that these downflows con- 
verge into downdrafts, forming tree-like structures down 
into the convection zone, where the intergranular lanes 
combine into more and more widely separated branches 
([Spruit et al.lfllMl 1 . 

As a result of upflows and downflows of plasma, the 
evolution of the cellular network in the photosphere 
is quite complex with granule s growing, decayin g and 
merging with other granules ([Spruit et al.l Il990f ). but 
some studies have been carried out trying to define and 
model some physical characteristics o f solar granulation 
like the size of granules (fo r example. iRoudier fc Mullen 
U987tlSchriiver et al.|[l997t ). lifetimes (ISpruit et al.lll990l) 
and flow profiles of the plasma (see ICanfieldlll976l and 
iDurrant et al.lll979l for estimates of vert ical velocity com- 
ponent. and iBeckers fc Morrison! 119701 for the horizon- 
tal velocity compone nt). The differences betwee n stellar 
convection zones (see lNordlund fc Dravinsll 19901 for some 
simulations regarding the convection in some solar-like 
stars), makes it likely that the structure of the granu- 
lation on other stars differs from the granulation pat- 
tern we have on the Sun. The change in magnitude 
and topology of the granular velocity field should affect 
the total Poynting flux transmitted into the corona, and 
most likely the amount of heat dissipated in stellar at- 
mospheres. 

Using a numerical code to simulate the granular mo- 

1 See page 291, Figure 7 for an illustration of this scenario. 
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tions at the photosphere, we follow the effects on the 
corona by measuring the electric current density. Run- 
ning the code for several values of parameters associated 
with the structure of the granular motions, we investi- 
gate what happens with the corona as a function of the 
granular characteristics. 

2. NUMERICAL PROCEDURE 

2.1. Method 

The simulations are carried o ut with the 3-D MHD 
code e xtensively described in iGudiksen fc Nordlundl 
(|2005bf ). which solves the fully compressible MHD equa- 
tions: 



dp _ 

-JL = -V.(puu- 5 ) 

Be 



VP + JxB + pg, 



— = -V • (eu) - PV • u 



Joule 



Qvi 



(i) 

(2) 
(3) 



dB 

1W = - VXE 



(4) 

E = ? ? J uxB, (5) 
/iJ = VxB. (6) 

Here p, u, P, J, B, g, /i, E and e denote the mass den- 
sity, velocity vector, gas pressure, electric current density, 
magnetic flux density, gravitational acceleration, the vac- 
uum permeability, electric field vector and the internal 
energy per unit volume respectively; The terms r and ry 
refer to the viscous stress tensor and the magnetic diffu- 
sivity; Q V i SC and Qj ou ie represent the viscous and Joule 
heating. 

The set of MHD equations is computed on a stag- 
gered mesh of 100 3 grid points corresponding to a box of 
20 x 20 x 20 Mm 3 (distance between two grid points of 
0.2 Mm), employing 6th-order derivative operators, 5th- 
order interpolation operators to align physical variables 
in space, and a 3rd-order Runge-Kutta time-stepping 
scheme. 

Simulations on this work use a classical heat con- 
duct ion along th e magnetic field, Spitzer conductiv- 
ity (jSpitzerl 119561 ) and also a cooling function repre- 
senting radiativ e loss es in the optically thin coronal 
plasma. iParkerl (|1972| ) argued that the vigorous forma- 
tion of tangential discontinuities (current sheets), mag- 
netic dissipation is practically independent of the con- 
ductivity. Moreover, i t wa s showed n umerically by 
iGalsgaard fc Nordlundl (|1996f ) and also iHendrix et al.l 
( 19961 ) that the energy created through dissipation is not 
strongly dependent on small scales. One should expect 
that, by resolving the smallest scales in a simulation of 
DC heating, the dissipated energy could have a different 
value. Instead, high resolution does not give dramatic 
changes in results as the stress applied on the magnetic 
field is not increased but rather shifted to a fine scale. 
The important point on the ohmic/nanoflare heating hy- 
pothesis is the level of stress experienced by the magnetic 
field which does not depend on the resolution scale, so 
consistent results can be retrieved from simulations with- 
out the necessity of using very high resolution. 

The MHD equations are solved in a 3D x-y-z co- 
ordinate system with periodic boundary conditions in 



the horizontal direction. The top of the box has zero- 
gradient of temperature at the boundary so no heat ex- 
changes are allowed. Also the vertical gradient of the 
horizontal velocity is null at the upper boundary. Both 
the density and potential field are extrapolated in the 
ghost zones (regions outside the computational domain) . 
Vertical velocity is set to zero on both top and bottom 
boundaries and the driver establishes the horizontal ve- 
locity at the bottom. The lower boundary is stressed 
by a time dependent velocity field, constructed from a 
Voronoi tesselation, reproducing the granulation pattern. 
The velocity field reproduces the geometrical pattern as 
well a the amplitude power spectrum of the velocity and 
the vorticity, leaving no free parameters. 

The photospheric motio ns are produced through a 
Voronoi tessellation driver (jOkabe et a,l.ll2000l ) which re- 
produce the geometry and flow pattern of the granular 
surface of the Sun. A pa rticular model o f a Vo ronoi tes- 
selation investigated by ISchriiver et alJ (|1997| ). the so- 
called "multiplicatively weighted Voronoi" tessellation, 
was found to reproduce the granular pattern of the Sun 
very well. It segments the two-dimensional space into 
Voronoi regions V(p i ) defined by 
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for j e [1, 



(7) 



with Wi representing the different weights (greater than 
zero) associated with the generator points x i: and p de- 
noting the index of the power law function which controls 
the tessellation. By using these weighted distances be- 
tween the generated points, the generation of very large 
or very small granules can be avoided. By choosing the 
weights Wi carefully, a uniform granular structure is pro- 
duced with only acceptable sizes for the granules, ranging 
in size between 1 and 2 Mm, in the case of solar photo- 
spheric granules , to ~ 30 Mm or larger for supergranules 
(|Leightodlf963l ). 

The velocity p rofile used by the model is taken from 
iSimon fc Weiss! ([1989), who fitted observational data 
from the SOUP instrument on Spacelab 2, and then ad- 
justed to better fit the solar velocity spectrum by 



Vh(r) = v,(t)r 2 e \ r °< 



(8) 



In the previous equation Vi is a time-dependent quantity, 
proportional to Wi (t) which is defined by 



Wi(t) = Vt^exp 



t - tr 



(9) 



with Wi being the maximum weight and q an even integer 
related with the growth time of the granule. The velocity 
at the center of granules evolve in time along a Gaussian 
profile Vi cx e~(*/ to )". 

2.2. Simulations 

To quantify the effect of the photospheric granular net- 
work on the corona, we set up three sets of simulations. 
In each set we vary one of three parameters describing 
the granules: The granular radius r, the granular lifetime 
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TABLE 1 

Main parameters for the simulations 



Input parameters 

Simulation Set r (Mm) r (min) v (km ■ s _1 ) # Snapshots t (min) 



I 1, 2, 3 2 5 

II 1 3 5 

III 1 5 5 

IV 2 2 8 

V 2 2 12 

VI 3 2 5 

VII 3 2 5 



3 600 259.6 

3 600 251.3 

3 300 128.7 

3 750 318.7 

3 1000 414.1 

6 600 261.3 

10 750 326.8 



Note. — Variables on the table are r, the size of the granules; r, the lifetime of the 
granules; vq, the plasma flow velocity, t, the time-span of the simualtion. 



r and the velocity amplitude of the plasma Vq. Each 
set contains three simulations with one simulation being 
identical for all the sets so the total number of simula- 
tions is 7. The details of the simulations and sets are 
provided in Tabic Q] 

Simulation I is common to all the three sets and cor- 
responds to the closest reproduction of solar conditions, 
by choosing r = 2 Mm, r = 5 min and vq = 3 km s _1 . It 
is important here to point out that these parameters are 
not entirely independent and setting these values does 
not produce granules that have exactly these character- 
istics. The actual values for the granular radius, granular 
lifetime and average velocity in the granular field will be 
discussed in Sec. 14.11 

In set 1 both the lifetime and velocity amplitude are 
kept constant, while the granular radius is varied from 2 
to 5 Mm. In the second set we simulate granules with 
lifetimes of 5, 8 and 12 minutes, maintaining the size of 
the granules and velocity amplitude constant and equal 
to 2 Mm and 3 km s _1 respectively. The third set in- 
cludes simulations with three different flow velocities, 3, 
6 and 10 km s _1 and granule sizes and lifetimes set to 2 
Mm and 5 minutes, respectively. 

Our simulation domain is initially a current-free space 
containing vertical magnetic field lines extending from 
the base which coincides with the photosphere (z = 0) 
to the top of the box, 20 Mm up in the corona (z = H = 
20 Mm). We then let the code run for several hundreds 
of minutes. The temporal resolution, i.e. the interval be- 
tween two consecutive snapshots produced by the model 
is roughly half a second, thus providing us with an al- 
most continuous evolution of the desired variables. We 
should note that the time span of each simulation is de- 
pendent on how fast the electric current density reaches 
a steady-state. After this point (designated here as the 
turn-over point) we extend the simulation, in order to 
ensure that the steady-state was reached. 

3. RESULTS 

We are only interested in the behaviour of the dis- 
sipated current in the corona. For all the simulations 
we follow the average current density in the corona be- 
cause it is a direct measure of the total dissipated mag- 
netic energy. To avoid effects of the timevarying local 
height of the transition region, we calculate the aver- 
age current density, taking into consideration the vol- 
ume between z = 5 Mm and z = 15 Mm. The size 
of our box is small compared with the typical size of 
coronal loops. The purpose is to study a portion of the 



loop from just above the transition region up to several 
mega-meters into the corona. This will correspond to 
one of the loop legs where the heating becomes more im- 
portant. The magnetic field becomes nearly force- free 
above the transition region and the current along each 
field line will b e proportional to the magnetic field and 
as described bv lGudiksen fc Nordlundl (|2005a( ) the aver- 
age heating (which is proportional to the magnetic field 
strength squared) drops exponentially above that layer. 

One behavior is common to all the simulations. The 
evolution of the average current density in the corona 
shows a very steep increase at the beginning of each sim- 
ulation, followed by a gradual deceleration and a turn- 
over point, and later a stabilization around a constant 
level designated here as asymptotic average current den- 
sity squared (J\). Nonetheless, some differences exist 
between the experiments. 

To better identify and understand these differences we 
perform a curve fit for the data assuming a special vari- 
ation of the logistic function, given by 

F(t) = a -e-^ (10) 

where a tl i = 0, 1 are the two free fit model parameters. 
One important property of this function is that it has a 
horizontal asymptote, i.e. F asyrn = lim F(t) = ao, so it 

t— *oc 

can be used to determine the asymptotic level of the aver- 
age current density when the box reaches a steady-state. 
The turn-over point (t turn ) is defined here as the time 
when the fitting-function reaches 80% of the asymptotic 
value of the fitting-curve: 

Fturn — 0.8 ■ Fasym (H) 

The turn-over points as well as the J\ and the fitting 
parameters are summarized in Table [2j In particular, re- 
sults from the simulations are plotted in Figures Q] - 31 
The figures clearly illustrat the fast initial growth of the 
average current density, strongly impelled by the sudden 
advection of the magnetic field lines. The random char- 
acter of the movements of the footpoints is the natural 
way to interpret the fluctuations of the average current 
density even after the turn-over point. 

Even though the granular pattern is changing, the sta- 
tistical properties of the velocity field are constant, so 
the sharp rise in the average current must mean that 
the magnetic field is initially able to contain the energy 
injected into it by the Poynting flux produced by the 
photospheric velocity field. At the turn-over point a sat- 
uration is reached, and the magnetic field dissipates as 
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TABLE 2 

Fitting parameters (an. and ai) in Eq. [To] their 

UNCERTAINTIES (e ao AND £ ai ) AND TURN-OVER POINTS. 

Fitting-function 



Simulation 


Turn-over point (min) 


a 


£a 


ai 


e ai 


I 


54.5 


6.23 


0.06 


12.2 


0.7 


11 


36.7 


7.93 


0.08 


8.15 


0.61 


III 


25.7 


11.6 


0.2 


5.78 


0.49 


IV 


95.4 


8.76 


0.11 


21.2 


1.3 


V 


150.0 


10.97 


0.15 


33.2 


2.0 


VI 


61.5 


6.80 


0.06 


13.8 


0.7 


VII 


81.5 


7.50 


0.08 


18.2 


1.0 



much energy as is injected into it by the Poynting flux 
at the photosphere. 

At first sight, the evolution of the average current den- 
sity in the corona is quite similar for all the experiments, 
despite the variations of r, t and vq. Thus, if we do 
not consider the numbers involved but concentrate only 
on the behavior of the average current density, we con- 
clude that it is practically impossible to distinguish be- 
tween different simulations. After the exponential in- 
crease of the average current density, the magnetic field 
line mapping inside the box starts to become quite com- 
plicated and the shear pattern will evolve independently 
of the parameters which govern the motions of photo- 
sp heric material. A similar arg umentation is presented 
in lGalsgaard fc Nordlundl (|1996f ) with respect to the vari- 
ation of the driving boundary pattern. 
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Fig. 1. — Evolution of the average current density in the box as 
a function of time (solid line) and its curve fitting function (dash 
dotted line) in simulation I with granules size/ granules lifetime/ 
flow velocity = 2 Mm/ 5 min/ 3 km s _1 . Also shown are the J 2 A 
(horizontal dotted line) and the turn-over point marked with the 
small filled circle over the fitting curve. The bottom part of the 
plot shows the residuals in the sense of (F-C: Fitted-Computcd). 



4. DISCUSSION 

To better perceive the consequences of different photo- 
spheric configurations on the average current density, we 
analyze the moment of turn-over point and the J\ at the 
steady-state as wc change the three variables: granular 
size, granular lifetime and flow velocity of the plasma. 
This is done by having in the background the theory 
proposed by Parker for the coronal heating, based on the 
dissipation o f energ y through electric currents (|Parkerl 
11973. ll983al lH ll988l ). Assuming a similar box with ver- 
tical magnetic field at t = 0, he argued that the photo- 



spheric motions due to granulation drag the footpoints 
of the magnetic field, stressing it and creating current 
sheets which will dissipate and heat the corona. 

Taking the velocity of the photosphcric motions (equiv- 
alent to the shuffling velocity of the magnetic field lines) 
as v, the angle of deviation of the magnetic field as 8, 
time as t and vertical size of the box represented by H, 
the following equation, 



(12) 



is valid for relatively small values of 8. By inserting the 
the vertical component of the magnetic field B, and the 
transverse component, B± = B tan 9(t) we find that 



Bj 



Bvt 
If' 



(13) 



The onward motion of the footpoint creates tension in 
the field line in the opposite direction, so the movement 
of the footpoint does work on the magnetic field and the 
energy input grows as the field is dragged. After some 
time when 8 reaches a "critical angle", the dissipation 
of the stressed field exactly counters the produ ction of 
stress by the continuing motion of the footpoints. iParkerl 
estimated the critical angl e based on obse rvation and en- 
ergy arguments to be 20° (|Parkerlll983at ) or 14° (|Parkerl 
[19881) . 

4.1. The driving velocity field 

The details of the driving velocity field in the photo- 
sphere is important. As seen, all the three parameters 
controlling the velocity field have an impact on the J\ 
and the turn over point. 

In a simplified picture, the granular velocity field will 
advect magnetic footpoints radially in each granule, un- 
til they arrive at the intergranular boundary. Depending 
on the location of the surrounding granules, the mag- 
netic footpoints will then cither stay at the location in 
the intergranular lane or be advected almost perpendicu- 
lar to their previous velocity until they reach an intersec- 
tion between more than two granules, where the velocity 
is generally low. The development of new granules will 
make that process repeat itself. The direction of move- 
ment will be independent of the previous history of the 
magnetic footpoints, so they will experience a random 
walk, where the step length will be roughly the granular 
radius, unless the velocity of the flow inside the granules 
is so low that they will never have time to reach one gran- 
ular radius before the granule disappears. In that case 
the random walk step length will be the product of the 
velocity inside the granule and the lifetime of the gran- 
ule. A convincing argument that the magnetic footpoints 
do experience a random walk, is found by following mag- 
netic footpoints in the granular velocity field, assuming 
they are passively advected (see Fig. [5] for an example of 
the path travelled by one footpoint). We have followed 
1024 footpoints at the surface (forming a 32 x 32 grid at 
the bottom of the box) over time for each simulation. We 
then plot the displacement from the original position (R) 
as a function of time. Figure [6] shows that total displace- 
ment of the footpoints in the granular field of simulation 
I, and a fit of the form 



R = Ct p . 



(14) 



50 100 150 200 250 20 40 60 80 100 120 
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Fig. 2. — Evolution of the average current density in the box as a function of time (solid line) and its curve fitting function (dash dotted 
line), in simulations II (left) and III (right). 
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Fig. 4. — Evolution of the average current density in the box as a 
line), in simulations VI (left) and VII (right). 

In the previous expression C and p are constants and a 
value of p = 0.5 would correspond to theoretical random 
walk. The fitting curves have values of p as indicated on 
Table[3]which are fairly close to the theoretical value. We 
will assume that the motion experienced by a magnetic 
footpoint is a random walk, and we can then fit the data 
with 

R = eVN. (15) 

The random walk is described by two parameters: the 
number of steps N and the step length I. The two pa- 
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:tion of time (solid line) and its curve fitting function (dash dotted 

rameters cannot be found independently, so they are not 
easy to find as we have only snapshots in time of the sim- 
ulation. They can be found indirectly if we can find the 
step length I and the step time f (the time required for 
one step) . If these can be found the number of steps can 
also be determined since N = t/f where t is the elapsed 
time. Because we have decoupled the velocity inside the 
granules, their lifetimes and their size, we have to be 
careful with the estimation of the parameters. The step 
length is assumed to be the smallest of these two values: 
the product of velocity times lifetime or the average ra- 
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TABLE 3 

Values of p and 

UNCERTAINTY t v IN EQ. [14] 

FOR EACH SIMULATION 



Simulation p e p 



I 


0.630 


0.007 


11 


0.311 


0.005 


111 


0.610 


0.013 


IV 


0.547 


0.006 


V 


0.530 


0.005 


VI 


0.585 


0.009 


VII 


0.623 


0.005 



dius of the granules at the surface, i.e., the distance from 
the center of each granule to the intergranular lanes/ 
vertices. By making this assumption we assume that the 
footpoint starts its movement always at the center of a 
granular cell which is a strong simplification and it will 
be discussed later. A footpoint can initiate the motion 
anywhere inside the granule migrating then towards the 
intergranular lanes, probably lying there for a while until 
a new granule emerges and propels it again. So, the gran- 
ular radius should be understood as an over-estimated 
step length and the maximum value that this parameter 
can take. 

The average radius of the granules r av is 

calculated using the same approach taken in 

iWedemever-Bohm fc Rouppe van der Voord (|2009f ) 
for solar observations, in this case from a power spec- 
trum of the velocity field on the z — plane (see Fig. 
[7] for an example of the spectrum as a function of 
wavenumber for one snapshot of simulation I), where 
the value on the x-axis corresponding to the peak of 
the curve provides the diameter of the granules. From 
Fig. [7] it is immediately seen that the power spectrum 
has a broad peak, around D ~ 2.8 Mm, which gives 
r av ~ 1.4 Mm for this simulation. 

The average velocity v inside the granules is then cal- 
culated from the velocity profile (Eq. [8] and Fig. [8]) using 
the average radius from the powerspectrum of the veloc- 
ity field: 

l 

v = / Vh dr (16) 

r av JO 

The average lifetime f of the granules is found by look- 
ing at the distribution of the number of granules accord- 
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Fig. 6. — Distance R from the original location as a function 
of time for 1024 footpoints in simulation I. The dash dotted line 
corresponds to the average distance and the solid line represents 
the fit assuming a fitting function R = Ct p . 

ing to their development. Figure [9] shows an example 
(simulation I) which shows that most of the granules are 
living within 3.5 minutes around the instant of maximum 
development (t—toi), which will correspond to an average 
lifetime of 7 minutes. 

Initially £ is set equal to the product vf. In case the 
product vf is greater than r av , r av is taken as the step 
length. 

Table [4] organizes all these quantities discussed above 
for our simulations. 
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Fig. 7. — Velocity power spectrum vs. wavenumber for a snap- 
shot of simulation I. The figure shows the peak as a sign of the 
granulation pattern at the bottom of the box. 
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Fig. 5. — Random motion of the footpoint of a magnetic field 
line, at the bottom of the box. The square symbol marks the initial 
position and the 'X' shows the position of the footpoint at the end 
of the simulation I. 



4.2. Turn over point and critical angle 

Figure [10] (right) shows the correlation between the 
granular size and the time of the turn-over point of the 
simulations. A plausible explanation of this behavior 
consistent with Parker's model and the random motion 
of the footpoints relies on the critical angle and on the 
fact that footpoints must be moved a certain distance 
from their original location in order to reach that an- 
gle. The properties of the random motion means that 
the total distance traversed by the footpoints is more 
sensitive to the step length than the number of steps. 
All simulations of set 1 (where the size of the granules is 
changed keeping the other initial parameters constant) 
have approximately the same step length but one has to 
be careful as the size of granules affects the results. We 
have been assuming the center of the granule as the orig- 
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Fig. 8. — Example of the velocity profile of a granule belonging 
to simulation I. The velocity can never be zero inside the code in 
order to keep the photospheric motions, so instead of the dashed 
portion of the plot a constant value of v = 0.2 km s _1 is assumed. 



Fig. 9. — Histogram showing the distribution of granules ac- 
cording to their development in simulation I. Negative values on 
the x-axis correspond to the period pre-full development while the 
positive values represent post-full development phase. 



TABLE 4 

Average granular size (r a „), step length (£), average 

VELOCITY (v), AVERAGE LIFETIME (f) AND AVERAGE NUMBER OF 
RANDOM STEPS (N) FOR EACH EXPERIMENT 



Simulation r av (Mm) I (Mm) v (km ■ s 1 ) r (min) N 



1 


1.26 


1.10 


2.63 


7 


37.1 


11 


1.66 


1.10 


2.29 


8 


31.4 


III 


2.63 


1.11 


2.17 


8.5 


15.1 


IV 


1.32 


1.32 


2.58 


11 


29.0 


V 


1.29 


1.29 


2.63 


17 


24.4 


VI 


1.25 


1.25 


5.13 


7 


37.3 


VII 


1.30 


1.30 


8.61 


7 


46.7 



inal location of footpoints which is not necessarily true of 
course. Eventually a footpoint may start to be dragged 
very close to the intergranular lane so it cannot travel the 
distance equivalent to the step length. In this first set of 
simulations I ~ 1.10 Mm which means that for a foot- 
point to be able to move this distance it has to be located 
at least one step length from the intergranular lane. If 
the footpoint starts to move from a point farther away 
from the intergranular lane, its motion will be limited to 
I by the lifetime of the granule. There will be then an 
area of the granule around the center where footpoints 
can start the motion with high probability of traveling 
along £. This area inside granules is obviously larger and 
larger as we increase the size of the granules meaning 
that, statistically, the real step length will closely match 
the value found for larger granules and will progressively 
reduces as we decrease the granular size. Our values 
of step length correspond in this way to the maximum 
value the step length can take. Taking the size of the 
granules into account it is reasonable to think then that 
our value of the step length will diverge more and more 
from the real value as we decrease the size of the gran- 
ules. This can justify the inverse correlation observed 
between the granular size and the time for the turn-over 
point: smaller granules will have a smaller and smaller 
value of step length and then the equilibrium between 
current generated and released in the corona will happen 
at a later time when compared with larger granules. 

Figure [TT1 (right) shows the variation of the instant of 
turn-over point as a function of granular lifetime. Here 
the correlation is linearly increasing. According to the 



random walk theory there should be no dependence on 
the granular lifetime. Nevertheless we should not ex- 
clude the influence of the adjacent granules and their 
lifetime. For shorter lifetimes the topology of velocity 
field changes faster, and the magnetic field lines trapped 
in intergranular lanes will be moved sooner due to the 
faster changing velocity field topology. A second effect is 
that the electric current density is lower as we will see in 
Sec. [O] 

The time of turn-over point is also influenced by the 
velocity field and the correlation is shown in Figure [12] 
(right). It is obvious that as we increase the parameter 
vq in the simulations, the average step velocity (v) also 
increases approximately by the same amount, which will 
induce an increase of the step length. Associated with 
this fact, the electric current produced also dramatically 
increases and the time needed to balance this process by 
dissipating the larger amount of current will be longer 
than for low velocities. 

At this point there are enough elements to give an es- 
timate (in a rough way as we are doing some strong as- 
sumptions) for the critical angle introduced by Parker 
in his model. This critical angle should correspond to 
the average deviation from the vertical of the magnetic 
footpoints when the dissipation starts to cancel the trans- 
verse component of the magnetic field at the same rate it 
is created by the drag of footpoints in the photosphere. 
In this study that moment will coincide approximately 
with the time of turn-over point, when after the fast ini- 
tial growth of electric current density we observe a stabi- 
lization of this quantity caused by the balance between 
creation and dissipation of electric current. 

Using the values of the time of turn-over point, time 
step and step length and applying them on Eq. [15\ one 
can calculate the average displacement of the magnetic 
footpoints from its original location 

Rturn = i^^ (17) 

where t turn is the time of turn-over point and f the time 

Step, SO tturn/ f = N turn . 

The estimated deviation angle from the vertical will 
then be 

9 est = arctan ( I (18) 
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with H being the height of the box. H = 20 Mm. Results 
are indicated on Table 

Due to all the assumptions taken along this work, it 
is hard to find correlations between the angle and the 
initial parameters, or at least to put a high level of trust 
in them. Nevertheless one can observe that there are no 
huge variations of the estimated angle which varies from 
5.5 to approximately 12.5, thus leading to an average 
value around 10 degrees. This is not very far from the 14° 
postulated by Parker and could mean that our assump- 
tions are justifiable. There is a particular discrepancy 
in set 1 of experiments where the estimated angle de- 
creases from 8.76° in simulation I to 5.50° in simulation 
III. According to fTS] the critical angle is proportional 
to the displacement at the moment of turn-over point, 
which is proportional to the step length from [171 We 
should expect a smaller angle for smaller granular size as 
the turn-over point is reached faster leading to a shorter 
displacement and the calculated step length is about the 
same in all these 3 simulations. However it was already 
explained at the beginning of this section that the true 
step length should be smaller for smaller granules, so this 
fact would uniform the critical angle, making it approx- 
imately the same for the three simulations. 

4.3. Average electric current in the atmosphere 

The energy dissipated by currents in the atmosphere 
originates as Poynting flux injected into the magnetic 
field by the photosphcric motions. The amount of dissi- 
pated energy in the atmosphere must therefore be pro- 
portional with the injected energy flux when a statistical 
equilibrium is reached. If this is correct there should ex- 
ist simple relations between the granular paramet ers and 
the ave rage electric current in the atmosphere. iParkerl 
(|1983af ) estimated the Poynting flux injected into the at- 
mosphere: 

vB± B_ BWAt 
4tt AttH v ' 

where v is the velocity which the magnetic footpoint is 
moved with (measured in km/s) during the time interval 
At, B is magnetic field strength, B± is the field strength 
of the transverse component of the magnetic field gen- 
erated by the movement of the footpoint and H is the 
height of box. As there is only a dependency on the flow 
velocity, the relations should here be rather simple. It 
is though important to notice that the velocity here is 
the effective velocity, meaning the velocity that carries a 
magnetic footpoint further away from it initial position 
in a straight line. It is therefore not that simple to es- 
timate the effect of the variables defined earlier, r av , f, 
and v. 

Figure [TU] (left) shows the dependence of the average 
current density in the atmosphere as a function of the av- 
erage granular size and the good linear fit between these 
two quantities: the current density generated goes up 
as the average size of the granules increases. The step 
length calculated for set 1 is roughly the same but be- 
cause of the reason already pointed out in the previous 
discussion about the variation of the instant of turn-over 
point with granular size, one should expect a smaller 
value of the step length than the calculated value for 
smaller granules. So in a photosphere with larger gran- 
ules the motions of the plasma will cause a higher dis- 



TABLE 5 

Values of the average 
displacement at the turn-over 

POINT (Rturn) AND ESTIMATES FOR 
THE CRITICAL ANGLE(0 es £ ) 



Simulation 


Rturn (Mm) 


dest 


I 


3.08 


8.76° 


11 


2.35 


6.71° 


111 


1.92 


5.50° 


IV 


3.89 


11.0° 


V 


3.83 


10.9° 


VI 


3.71 


10.5° 


VII 


4.44 


12.5° 



placement in a straight line of the footpoints of the mag- 
netic field lines, twisting the magnetic field more than 
in the case of smaller granules and thus producing more 
electric current. With the module of the Poynting vector 
proportional to the transverse component of the mag- 
netic field, B±, a larger granule, produces higher B± by 
creating a stronger twist of the magnetic field and there- 
fore produces an increase of the energy flux. 

Figure fTTI (left) evinces a correlation between the aver- 
age granular lifetime and J\. Magnetic footpoints caught 
inside a granule that lives longer follow a straight path ra- 
dially inside the granule for a longer time, and therefore, 
the effective velocity is higher. This happens as long as 
the granules arc large enough for the magnetic footpoints 
not to reach the edges of the granules within the granular 
lifetime, but in this set 2 of simulations the granular size 
is kept constant and then the influence of the size of the 
granules is the same for the three experiments of this set. 
On the other hand, inside short-lived granules footpoints 
will move straight in shorter time steps decreasing the 
effective velocity and consequently the J\. 

In Figure[T2lfleft) we show an identical scenario for the 
variation of the J\ as functions of the average velocity 
field at the photosphere. From equation [191 it seems that 
there should be a correlation between the total current 
squared and the velocity squared, but again the velocity 
in equation [19] is the straight line velocity while the ve- 
locity experienced by the magnetic footpoints are not the 
effective straight line velocity. The correlation between 
the J\ and the velocity can then be explained by the 
random walk nature of the velocity field. 

One can think, based on expression [12] that a similar 
correlation should be obtained between the average cur- 
rent density or the turn-over point and the product vf. 
Following the previous discussed relation of average cur- 
rent density and the angle 9, the latter should increase 
with vf, and the same should apply for the turn-over 
point. Figure [TBI shows this behavior, although it is not 
possible to linearly fit all the data points. This non-linear 
increase of both the average current density and turn- 
over point with vf implies that these variables are not 
connected in a simple way, and that the relation between 
the amount of electric current arising in the corona and 
the deviation angle 9(t) to the vertical should be more 
complex. This is probably again an effect of the differ- 
ence between the v and the velocity used in equation ll9l 

5. CONCLUSION 
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(left) and turn-over points (right) as a function of the average granular size, from simulations of set 1. 
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Fig. 11. — (left) and turn-over points (right) as a function of the average granular lifetime, from simulations of set 2. 




5.5 [ 



y=a+bx ~ 

a=5.68213 " 
b=0.213415 . 



4 6 

Average flow velocity (km s" 1 ) 



10 



^ 60;- 



30 E 



y=a+bx- 

a=40.7726 - 
b=4. 59269 = 



4 6 

Average flow velocity (km s~ 1 ) 



10 



Fig. 12. — (left) and turn-over points (right) as a function of the average flow velocity, from simulations of set 3. 



At the present time, we have the opportunity to de- 
ploy significant computing power available to explore in 
a deeper way the heating process of stellar atmospheres. 
Because observations from stars with the same resolution 
we get for the Sun are impossible to obtain due to the 
distance, numerical simulations stands as the only way to 
explain some of the more complex behaviours observed 
on stars. 

In this work we perform a number of simulations 
with the intention of studying how photosphcric param- 
eters affect the amount of electric current density in the 
corona. In particular we analyze three parameters, the 



size of granules, lifetime and the velocity of the plasma 
at the surface, using a 20 3 Mm 3 box with a stellar quiet 
photosphere at the bottom. We have seen that tuning 
these three variables at the surface and thus, simulating 
different environments at the stellar photosphere, we get 
an initial sharp increase of the electric current density 
induced by the motions of the plasma dragging, twisting 
and shuffling the magnetic field lines. Later, the produc- 
tion of electric current decelerates, reaching a state with 
average constant value after the turn-over point. Oscilla- 
tions around the J\ are a result of the random behavior 
of the advection of the footpoints of the magnetic field. 
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Fig. 13. — (left) and turn-over points (right) as a function of 
I, IV, V, VI and VII. 

According to our simulations, the time when the elec- 
tric current density gets to the steady-state is a conse- 
quence of the balance between creation and destruction 
of B± as well as the value of the j\. Both of them 
changes as we change the parameters r, r and v. Based 
on our experiments, we conclude that an increase of any 
one of these parameters will lead to a higher value of 
the J\, and that this relation of direct proportionality 
can be described with a linear function. Higher values of 
electric current density at the steady-state are obtained 
for larger granules as a result of creation of B± because 
the linear radial velocity inside the granules can work 
longer due to the smaller risk of the magnetic footpoint 
hitting an intergranular lane. The step length, the dis- 
tance travelled by a magnetic footpoint every random 
step, plays an important role here as the production and 
release of current in the atmosphere is closely associated 
with it. The lifetime of granules is also a factor which 
affect the electric current density on the corona. Follow- 
ing Parker's idea, the deviation angle 9 is proportional 
to the time travelled by a footpoint along the same di- 
rection. Thus, longer lifetimes mean higher horizontal 
displacements of the magnetic footpoints which cause an 
increase of the asymptotic electric current density. The 
turn-over point is also affected by the lifetime of the gran- 
ular cells at the photosphere. With short-lived granules 
the surface of the star changes quickly, the footpoints of 
the magnetic field arc shuffled, twisted and braided faster 
culminating in earlier turn-over points. 

Parker's argumentation can again be used to explain 
results of the variation of flow velocity. We have obtained 
an increase of electric current density at the steady-state 
as the velocity goes up. As of [131 the transverse com- 
ponent of the magnetic field is correlated with the flow 
velocity, so higher flow velocities generate more electric 
current. The resulting Poynting flux is proportional to 
the driving velocity and at the steady state, that Poynt- 
ing flux must lead to a total dissipation that is identical 
to the Poynting flux injected into the magnetic field. 

By doing these experiments, and varying only three pa- 
rameters of the photosphere, we have shown that these 
parameters can affect the amount of electric current den- 
sity in the corona and consequently the energy dissi- 
pated. The correlations obtained show that the size of 
the granules on the stellar surface, their lifetimes and the 
magnitude of the velocity of the plasma in the granules 
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affect the amount of current in the corona, which have in- 
fluence on the process of coronal heating. In the context 
of electric current heating fed by the photosphere these 
parameters seem to assume particular importance by de- 
termining the way the magnetic field is tangled above 
the photosphere and how energy will be released. Our 
simulations allow us to make an estimate of the so-called 
critical angle, first introduced by Parker, to be around 
10°. Although some strong assumptions are made in this 
work, this va lue does not l ook very distant from the 14° 
sugg ested bv lParkeri (fl98l . 

However, we should also point out some limita- 
tions of our model. Although this model is consis- 
tent with the ideas supporting the nanoflare coronal 
heating hypothesis, we can not resolve single events. 
iGalsgaard fc Nordlundl (|1996l ) shows and argues that the 
resolution is not the important factor as long as it is rea- 
sonable, since the Poynting flux is the provider of the en- 
ergy, and as long as the injection mechanism is resolved, 
in this case the granules, the dissipation locations will 
be correct, the total dissipated energy will be correct, 
but the fine scale structure of the dissipation will not 
be resolved. The Voronoi tessellation splits the bottom 
of the box into mosaics representing the granules which 
are governed by certain parameters. This is a simplified 
picture of what is a real surface of a star. The model 
used for stratification in the atmosphere is also simple 
and some details and contributions for the coronal heat- 
ing could have been missed as a consequence of that. In 
the interpretation of the results the way the step veloci- 
ties and step times are determined (averages in time and 
space) may lead to uncertain values, resulting then in no 
definite step length as well. In spite of the problems in 
evaluating the variables we use for the interpretation of 
the data, the relatively easy explanations of the results, 
make the values seem viable. 

Some studies have been done aiming to extract em- 
pirical relations from observations of cool stars between 
quantities such as X-ray flux density or magnetic flux 
density and some parameters associated with stellar 
properties like the Rq s sby n umber or color index. In 
particular iNoves et al.1 (|1984l ) suggested some relations 
between stellar activity level and a parameter similar to 
the Rossby number. They associated the Ca II H+K 
emission with a funtion g{B — V) which approximates the 
convective turn-over time (ic)- Also lVaughan fc Preston! 
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(| 19801 ) showed a correlation between Ca II H+K index 
and the color index, demonstrating a wide spectrum of 
activity level on solar-like stars. On other hand, a power 
law relationship between soft X-ray and the Ca II H+K 
fluxes, independent o f the stellar color, was found by 
iSchrijver et al.l (|1992j ). Flux in the X-ray band of the 
spectrum are mainly originated from the coronae of the 
stars, as a result of coronal heating mechanisms. By tak- 
ing this studies and relations and combining those with 
Parker's model for coronal heating which relies on the 
motions at the surface of the star to originate the energy 
release on the top of the atmosphere, we believe that 
stellar properties and ultimately photospheric parame- 
ters have some role to play in all this process of coronal 
energy losses. 

As the results show, distinct photospheric parameters 
result in distinct turn-over points and amount of electric 
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current density on the corona. This work shows that the 
characteristics of coronae can be influenced by the sur- 
face of stars which is directly connected with the mag- 
netic field above the photosphere. 
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